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Abstract. The paper addresses the two-point correlations of electromagnetic waves in 
general random, bi-anisotropic media whose constitutive tensors are complex Hermitian, 
positive- or negative-definite matrices. A simplified version of the two-frequency Wigner 
distribution (2f-WD) for polarized waves is introduced and the closed form Wigner-Moyal 
equation is derived from the Maxwell equations. In the weak-disorder regime with an arbi- 
trarily varying background the two-frequency radiative transfer (2f-RT) equations for the 
associated 2x2 coherence matrices are derived from the Wigner-Moyal equation by using 
the multiple scale expansion. In birefringent media, the coherence matrix becomes a scalar 
and the 2f-RT equations take the scalar form due to the absence of depolarization. A parax- 
ial approximation is developed for spatialy anisotropic media. Examples of isotropic, chiral, 
uniaxial and gyrotropic media are discussed. 



Consider the electromagnetic wave propagation in a random dielectric. Let u(x, t) = 
(D(x, t), B(x, t))''' be the displacement-magnetic-induction vector field. Then the mutual 
coherence function is given by |3] 



(1) (u(xi,ti)ut(x2,t2)> = / e^("^--i)V^-(-i+-^)/2(U(xi,u;i)Ut(x2,u;2)>M^^^2 



where (■) is ensemble averaging, t = (ti + t2)/2,r = ti — ^2 and U(x, cu) is the frequency 
component of u at frequency u. Throughout, all vectors are by default column vectors and 
t denotes Hermitian conjugation. 

Radiative transfer theory [U IH El [IHl [131 [HI [151 [IZj has been traditionally carried out in 
space-time with one time variable [ti = ^2)- The main goal of this paper is to derive equations 
for the quantity directly related to (U(xi, u;i)U^"(x2, c<j2)). In particular, we obtain the two- 
frequency radiative transfer (2f-RT) equations for polarized light in the weak-disorder regime 
with an arbitrary bianisotropic background. The 2f-RT equations then determine, via ([T]), 
the two-space-time correlations of the electromagnetic wave in random media. 

Our approach is set up for the most general linear local, lossless electromagnetic materials, 
in which each of the field vectors E and H is coupled tensorially to both D and B. Such 
materials have been the subject of considerable recent interest. One reason is that they can 
be created as metamaterials, i.e. composites of more conventional materials in which E is 
coupled to D alone and H is coupled to B alone. 
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1. Introduction 



The two-frequency approach has been previously pursued in terms of wavelength-rescaled 
two- frequency Wigner distribution (2f- WD) in the case of a uniform background . In the 
present work, we introduce an alternative version of 2f-WD and derive the corresponding 
2f-RT equations for the associated 2x2 coherence matrices in the case with an arbitrary 
background. We give several examples for which the scattering kernels can be computed 
explicitly. We show that birefringence naturally leads to decoupling of the polarization 
modes and the absence of depolarization in such media. As a result, the 2f-RT equations 
simplify to scalar equations. 

In Section [2] we formulate the problem in terms of the straightforwardly defined 2f-WD 
and derive the two- frequency Wigner- Moyal equation in Appendix A. In Section [3] we analyze 
the problem for high-frequency waves in an arbitrary background bianisotropic medium in 
the absence of random fluctuations. This is the geometrical optics regime. In Section HJ 
we consider the weak-disorder regime where, in addition to the arbitrary background, small 
random fluctuations are present on the scale of the wavelength. In Section [5] and [6] we employ 
the multiscale expansion technique to derive the radiative transfer equation from the two- 
frequency Wigner-Moyal equation for the weak-disorder regime. We also derive a paraxial 
approximation for the polarized light in a spatially anisotropic medium. In Section [7] we give 
several examples of isotropic, chiral, uniaxial and gyrotropic media for which the scattering 
kernels can be explicitly calculated. We conclude in Section [H] with a discussion of the flnal 
expression of the mutual coherence in terms of the solution of the 2f-RT equations. 

2. Maxwell equations and Wigner-Moyal equations 

In this paper, we consider the electromagnetic wave propagation in a heterogeneous, loss- 
less, bi-anisotropic dielectric medium. We assume that the scattering medium is free of 
charges and currents and start with the source-free Maxwell equations in the frequency lj 
domain 
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where K is, by the assumption of losslessness, a Hermitian matrix [12] 

(3) K = 

with the permittivity and permeability tensors K*^, K'^, and the magneto-electric tensor 
|16] . The Hermitian matrix K is assumed to be invertible. The present formulation encom- 
passes the acoustic, electromagnetic and elastic waves so that the 2f-RT theory developed 
here can be extended to these waves without major changes. We choose D, B as the primary 
flelds because they are transverse (divergence- free). 

In an isotropic dielectric, K*^ = el, K'^ = fil, = 0. In a biisotropic dielectric, as 
well as K^, K'^ are nonzero scalars. A reciprocal chiral medium is biisotropic with purely 
imaginary = ix- The appearance of nonzero arises from the so called magnetoelectric 
effect [11]. Crystals are often naturally anisotropic, and in some media (such as liquid 
crystals) it is possible to induce anisotropy by applying e.g. an external electric fleld. In 
crystal optics, K'^,K'^ are real, symmetric matrices and = [3]. In response to a 
magnetic fleld, some materials can have a dielectric tensor that is complex- Hermitian; this is 
the gyrotropic effect. A magnetoelectric, bi-anisotropic medium has a constitutive relation 

2 



(j3]) with complex Hermitian K^, and a complex matrix satisfying the Post constraint. 
It has been shown that a moving medium, even isotropic, must be treated as bi-anisotropic 

mm- 

Writing the total field U = (D, B) we introduce the two-frequency matrix-valued Wigner 
distribution 



(4) 



W(x,p;a;i,a;2) = ^ J e'^P^^Ui (x + ^)ut (x - ^)rfy 



where Ui and U2 are the total fields at frequencies Ui/i and 072/^ respectively. The parameter 
i is roughly the ratio of the wavelength to the distance of propagation. In the present setting, 
i <^1. Correspondingly, we will replace u in hj uj/i. The 2f-WD is clearly equivalent to 
the 2-point function U1U2 via the inverse Fourier transform. 
Notice the symmetry of the Wigner distribution matrix 

(5) W"^(x, p;u;i,u;2) = W(x, p;t^2,c^i)- 

In other words, the right hand side of (jlj) is invariant under the simultaneous transformations 
of Hermitian conjugation f and frequency exchange Ui ■^002- 

In what follows we shall omit writing the arguments of any fields unless necessary. 

We put the equation in the form of general symmetric hyperbolic system [T7] 

(6) - ijV + K,d,^ (K-^U) = 
where Rj are the symmetric matrix given by 
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The matrices iTj,j = 1, 2, 3 are related to the photon spin matrices [2J. For ease of notation, 
we set L = K~-^. 

Let u' = (ijj\—ijJ2)li and Uj = (co'i+ct;2)/2. The 2f-WD satisfies the Wigner-Moyal equations 



£q, 



(7) iu'W = -p,R, / e'''^^L(q)W(p-^)rfq-- / W(p + ^)L(q)e*'i^^rfqp,R 



iq 



+ -RA^ / e^'^"'L(q)W(p-^)rfq+-9.^. / W(p + ^)L(q)e^'i'^dqR 



£q. 



u— W 



£q 



iq. 



+ J e^'^"'L(q)W(p - -^)dq - -d,^ / W(p + ^)L(q)e^''"'rfqR,- 

where L is the Fourier transform (spectral density) of L 



Lfx) 



e^^^'iL(q)dq. 
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The derivation is given in Appendix A. For a Hermitian L we have 

L(p) = Lt(-p), Vp. 

Clearly eq. ([7]) is related to the time derivative of the mutual coherence ([T]) with respect 
to the central time t while eq. ([7]) is related to the time derivative with respect to the 
differential time r. Our will focus first on eq. ([7]) and comment on the constraint posed by 
(IH]) in the Conclusion. The full analysis of eq. ([8]) requires substantially different treatment 
and will be presented elsewhere. However, we will discuss the constraint imposed by the 
leading order terms of eq. ([8]) in the Conclusion and its implication on the two-spacetime 
correlation. 



V2W = -R,d,^ [LW] + -d,^ [WL] R,. 



3. Geometrical optics 

In this regime, we let £ <^ 1 implying a small ratio between the wavelength and the scale 
of background heterogeneity which is comparable to the distance of propagation. 

Let us first simplify eq. ([7]) by expanding the expression in the power of i and neglecting 
0(£)-terms. The first two terms on the right hand side of eq. ([7]) reduce to i^^Vo — Vi where 

Po(p)W = ipjRjLW - iWLpjRj 

while the last two terms on the right hand side become 

iR,a.^. [LW] + 1( 

We employ the regular expansion W = W + i'Wi + ■ ■ ■ and substitute it into the resulting 
equation. The leading-order equation 

(9) VoW = 

can be solved as follows [17j . 

For a positive (or negative) definite L, L^^^ is well-defined and under the transformation 
L^/^ the matrix pjRjh is transformed into the Hermitian matrix L^/^pjRjL^/^ which has 
a complete set of eigenvectors and eigenvalues {^"} C M. Let {d'^'"} be the associated 
eigenvectors in the original vector space where the index a keeps track of the multiplicity. 
Let the eigenvectors {d°"'"} be normalized such that d°"'"'''Ld'^''' = S^j^r^a^- is easy to check 
that {e'^'""'' : e*^'" = Ld*^'"} are the left eigenvectors of pjRjL and they are orthogonal to 
{d'^'''(p)} with respect to the standard scalar product: 

(10) e'^'^td-.C = 

Clearly, the eigenvalues as a function of the wavevector p define the dispersion relations. 
For general bianisotropic dielectric, it is easy to check that the zero eigenvalue = is 
always an eigenvalue with the associated left eigenvectors 

10,2 _ 



;il) e'^'^(p) ~ ( - 1 , e^'^ip, - 



Since L is invertible, it follows that the null space of pjRjL is spanned by {d'^'^ = Ke°'^, d° 
Ke°'2}. The relations 1^ and ([U]) imply that d^''^,r ^ 0, are transverse vectors in the 
sense that they are orthogonal to the wavevector p. 
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Throughout the Enghsh indices represent the spatial degrees of freedom while the Greek 
indices represent the polarization degrees of freedom. The Einstein summation convention 
and the Hermitian conjugation are used only on the English indices. 

Define 

(12) D"'"^(p,q) = d"'°(p)d"''^(q)t 

(13) E-'"^(p,q) = e'^'"(p)e'^'f(q)t. 

The null space of Vq is the linear span of {D'^'"''(p, p), Vr, a, C, p}, denoted by 9Jlp, for each 
p ^ with the scalar product Tr [HtLGL] , H, G G 9Jtp. 
Then the general solution to ([9]) can be expressed as 

(14) W = ^iy„-^D'^'<(p,p) 

where W^(- are generally complex- valued functions. The matrices W'^ = [W^^^] , free of the 
English indices, are referred to as the coherence matrices. 

The constraint that the electric displacement D and the magnetic induction B are both 
divergence-free yields 

(±V,±V) • W = 
which, in view of the definition (jlj), is equivalent to 

(15) (±p^±pt)w = o. 

Hence by f lTT]) e'^'-'^W = and by ffTOj) W° = where is the 2x2 coherence matrix 
associated with the non-propagating mode = 0. This implies, by ffTOl) . that W is a 
transverse field. 

The 0(1) equation is 

(16) PoWi = luj'W + PiW - V2W 

which is solvable if the right hand side is orthogonal to the null space, 9Jlp, of Vq. The 
solvability condition for f|T6l) then leads to the governing equation for the 2x2 coherence 
matrices: 

(17) iuj'W^ - V^VT ■ VpW^ + Vpfi^ ■ VxW^ - C^W^ - W^C^"^ = 
where the depolarization matrix = [C|q,] is given by 

Using ffTOl) we can cast in the explicitly skew-symmetric (in ^, a) form: 

The details of the calculation is given in Appendix B. Note that eq. (fT7|) is invariant under 
the simultaneous transformations of Hermitian conjugation and frequency exchange. 
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4. Weak-disorder regime 



Now we consider the weak coupling regime with the permittivity-permeabihty tensor K 
given by 



;i8) 



K-^(x) = Lo(x) \i + V~ev{j) 



< 1 



where the Hermitian matrix Lq represents the slowly varying background medium and v^V 
represents the medium fluctuations. The small parameter i describes the ratio of the scale 
of the medium fluctuation or the wavelength to the propagation distance or the variability 
scale of Lq . 

To preserve the Hermicity of L the matrix V must satisfy 



(19) 



VtLo = LoV. 



We shall assume below that Lq is either positive or negative deflnite. A negative-deflnite Lq 
gives rise to negative index of refraction [TSl HH]. A nondeflnite Lq gives rise to complex- 
valued refractive index and hence a lossy medium. To flx the idea, let us take Lq to be 
positive deflnite. Our method applies equally well to the negative deflnite case. 

We assume that V = [Vij] is a stationary (statistically homogeneous) random fleld with 
the spectral density tensors # = [$ij>rm], ^ = [^ijmn] such that 



(20) 
(21) 

which implies 



{Vij{:S.)Vmn{y)) 



(22) 
(23) 



^ijmnip)5{p - q) 
%mri(p)5(P + q) 



Here and below * denotes the complex conjugation. In the case of real-valued V, $ = 
The spectral density tensors have the basic symmetry 



(24) 
(25) 

Eq. f[TI?l) implies that 

(26) 
(27) 



ijmn \tr/ 



ijmn \ 



-^O.ij^mnjKp) 
Lo,ij^mnjl{p) 



^mnijip), 
^mraj(p), 



Lo,lj^mnji{p) 
^0,lj'^rnnji{p) 



z^'W = ^p,R, I e^'i^''Lo(q)W(p-^)dq-^ J W(p + ^)Lo(q)e^^'''ciqp,R, 

+ IrA, I e^'''''Lo(q)W(p - ^)rfq + ^d,^ J W(p + ^)Lo(q)e^'''''c?qR,- 
+ -^p,R, 1 e^'^^^i;!^) V(q - q')W(p - |)rfq'rfq 

+ ^R,a,, 1 e^''^^i;(|)V(q - q')W(p - |)dq'rfq 

(28) +^9,^, 1 W(p + |)Vt(q - q')i;(|)e^''^ Wc^qR, 

where x = x/£ is the fast spatial variable and (Vt)(q) = V"''(— q) the Fourier transform of 
As in the geometrical optics we approximate the first four terms on the right hand side 
of eq. ([28]) by i'^VoW - V{W + PaW. For the last four terms on the right side of eq. ([28]) 
we have 

y'e^q^*My)V(q-q')W(p-|)rfq'^;q ^ Lo(x) ^ e^'''*V(q)W(p - |)rfq 

j W(p + |)Vt(q - q')L;(^)e^^^*dq'ciq ^ ^ W(p + |)Vt(-q)e^''^*ciqLo(x). 
Hence we have the simplified form 

(29) iuj'w = r^VoW -ViW + r2'w + e-^/'^QiW + e-^/'^Q2'w 

where 

QiW = zp,R,Lo j e''i'*V(q)W(p-|)rfq-? J Wip + ^)\\q)e'^'^dqLoPjRj 

Q2W = ^R,Lo9,^y"e^'i'*V(q)W(p-|)rfq+^9,^Yw(p + |)Vt(-^^ 

Hereafter we shall work with eq. ( l29i) to derive the 2f-RT equations by emplying the multi- 
scale expansion (MSE) . 

5. MULTISCALE EXPANSION 

The key point of MSE is to separate the fast variable x from the slow variable x and make 
the substitution 

vw ^ v^w + rVxW. 

Consequently, 

P2W ^ PsW + r^PaW 
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with 

(30) 

(31) 



P2W 
P2W 



^R,9.^. [LW] + ia.^. [WL] R, 



The idea is that for sufficiently small ^ the two widely separated scales, represented by x 
and X respectively, become mathematically independent. 

We posit the expansion W = W + v^Wi + £W2 + substitute it into eq. ^ and 
equate terms of same order of magnitude. 

The 0(£~^) equation is 

(32) {V2 + Po)W = 0. 

We hypothesize that the leading order term W = W(x, p) be independent of the fast variable 
X. Thus 7^2 W = and eq. (1521) reduces to ^ and its solution takes the form f[T^ . 
The 0(£~"'^/^)-equation is 

(33) PoWi + P2W1 = -QiW - Q2W 

or equivalently, after Fourier-transforming- in-x and adding a regularizing 0(£)-term 
-i2£Wi(k, p) + A;jRjLoWi(k, p) + Wi(k, p)Lo%Rj 



+2 



PjRjLoWi(k, p) - Wi(k, p)LoPjRj 



(34) - 
and posit the solution 
(35) 



R,Lo(2p, + A;,)V(k)W(p - |) + W(p + |)Vt(-k)(2p, - A;,)LoR,- 



Wi(k, p) = J] C:^(k, p)D'^'"^ (p + ^, P - ^) 



a,a,C 



where C^^ are generally complex numbers. Note that the two arguments of D*^'"^ in (135|) are 
at different momenta p + k/2, p — k/2. 

We substitute ([H]) and ([35]) into eq. ([M]) and multiply it with e'^'"(p + 1)+ from the left 
and with e°"''' (p — |) from the right and solve the resulting equation algebraically. This yields 
the coefficients 

-1 



C:c(k,p) 



-n-(p + - |)e-°(p + |)tv(k)d-^(p - 1; 



(36) 



- \)w:^{v + |)d-np + |)tvt(-k)e-^(p - 1) 



When the leading term W is invariant under the simultaneous transformations of Hermitian 
conjugation f and frequency exchange uJi ^ U2, so is Wi which is equivalent to 

Ccai-K P; ^1, ^2) = C^(;(^, P; (^2, UJi). 



Finally the 0(l)-terms yields the equation after regularization 

(37) £W2 + ^R,Lo J-W2 + ^ J-WsLoR,- + tp.R.UW^ - ^WsLop.R, = F 
with 

(38) F = zcj'W + PiW-PsW- QiWi - Q2W1. 

It suffices to note that in order for the resulting solution to vanish in the limit £ — > 0, 
F must satisfy the solvability condition 

(39) limTr (G"fLoFLo) = 

for all random stationary matrices G satisfying eq. fl32|) . This can be seen by transforming eq. 
(I37D into Tr(GtLo(l37DLo) which by eq. ([32]) implies 2£Tr (GtLoW2Lo) = Tr(GtLoFLo) 
and hence fl39l) . 

Fortunately, we do not need to work with the full solvability condition fl39|) . It suffices to 
demand fl39l) to be fulfilled by all deterministic G, independent of x, such that 

(40) PjRjLoG - GLoPjRj = 0. 

In other words, as in (|T^ . we consider only a subspace of the null space of eq. fl32|) and 
replace fl39l) by 

(41) lhnTr(E-'«'^t(p^p)(F(x,i,p))) =0, Vr, e, J^, x, 5c, p 

where E'^''^'^ are defined in (fT3l) . As (1291) . (1331) and (1381) are invariant under the simultaneous 
transformations of Hermitian conjugation f and frequency exchange ui ■^002, therefore eq. 
(HTj) must also be invariant under the same transformations. 

To summarize, we have constructed the asymptotic solution W + -\/£Wi + £W2 which sat- 
isfies approximately the 2f Wigner-Moyal equation in the sense that the remainder vanishes 
in a suitable sense as £ — >• P, T]. 

With (!35l) - (l36l) and (!38l) . eq. (14T!) is an implicit form of 2f-RT equations that determines 
the leading order coherence matrix. Our next step is to write (14T!) explicitly in terms of 
explicit, physical quantities. 

6. 2F-RT EQUATIONS 

As in the geometrical optics, the terms icu'W + ViW — P2W in the expression (155]) yield 
the left hand side of (IT7]) after the operation (11T]) . 
First note the key expression 

I dqe^'itxv(q)Wi(p-|; 
J2 J dk (fi"(p + k) - n^ip) - ti)-^ 

X [ - n'^ip + k)iy;^(p)e7"*(p + k)vi/^,,,,(k)rf^-'''(p)D5"^(p + k, p) 

+Q'^{p)W^^{p + k)d';'^*{p + k)$*^,,(-k)e^-'^(p)D-'"^(p + k, p) 



and 



J (iqe^'itxw,(p + |)V(-q)t 

J2 I dk{n^{p)-n^p-k) 



S3 



X 



n'^ip)W;^ip - k)e7"(p)*<l>^,,„(k)rf^-'''(p - k)DZHp, p - k) 
+Q^ip - k)iy:^(p)d-'^*(p)vl/*^^„(-k)e7^(p - k)DZHp, P - k)" . 

The above expressions are independent of the fast variable x so (Q2W1) = 0. We also have 
TV(E-'«-t(p^p)(g^Wi)) 

= ^Y./ dk{n^{p + k)-n^{p)~tir'n'^{p)ef*{p) 

Q'-ip + k)W;MeT{p + k)^f,,,{k)d'^^'^{p)dr{p + k) 
+n'^{p)W:^{p + k)d'^^'^*{p + k)<|.}^,,(-k)e^''^(p)<'"(p + k)" 

dk{n'^{p)-n^p-k)-zi)-'n'^{p)e';^-'{p) 
X [ - i]-(p)iy;^(p - k)e7«*(p)<i.^,,„(k)d-'''(p - k)<'^*(p - k) 

+n'^(p - k)iyf^(p)t/-'^*(p)v[/*^^.„(-k)e7^(p - k)C^*(p - k)^ . 

To state the full result in a concise form, let us introduce the following quantities. Define 
the scattering kernel tensors ©"^(p, q) = [5^j^„^(p, q)] as 

(42) 5[„,^(p,q) = fi-(p)fi^(q)e:'«*(p)<'"(q)$..;,(p-q)e7'^(p)dP^*(q) 
Using fl2^ - (p71) one can derive the alternative expressions in terms of ^ 

5[„k(P' q) = ^^(P)^^(q)e?'^(p)^;'^*(q)*},.,(q - p)er (q)<'«*(p) 

= fi^(p)fi^(q)el'«*(p)<'"(q)vl/,,^,(q - p)e}'^*(q)rfp'^(p) 

and the properties 

(43) '5[,.c(P'q) = '^^^(P'q) 

(44) Sl^uM,p) = Sl,^^{p,q) 

For any 9Jtp-valued field G(p) define the z/)-component of the tensor ©"^(pjq) : G(q) 

as 

[6-(p,q) : G(q)]^, = 5[,,^(p, q)G„c(q). 

Define the tensors S"^ = [S^^^] analogous to the total scattering cross section as 

S-(p) = |y"5(fi-(p)-l]-(q))-z^ (n-(p)-r]-(q))-^|e-(p,q):Mq. 
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The 2f-RT equations for the coherence matrix then reads as 

(45) iuj'W^ - Vx^^^ ■ VpW^ + VpVT ■ VxW^ - C^W^ - W^C^"^ 

= 27ry"5(^]-(q)-^]-(p))6-(p,q):W-(q)rfq-S-(p)W^(p)-W^(p)S^t(p)^ Vr. 

The (5-function and Cauchy singular kernel arise because of the fact 

(46) lim = iTiSix) -\ — 

X -il X 

in the sense of generalized function. 

With the property (l43l) one can verify directly the invariance of (H5|) with respect to the 
simultaneous transformations of Hermitian conjugation and frequency exchange. 

6.1. Birefringence: scalar 2f-RT equation. Although, in view of flTT]) . the zero eigen- 
value n° = has multiplicity two in general, the nonzero eigenvalues in media other than 
the simplest isotropic medium often have multiplicity one as we shall see in Section [71 This 
is closely related to the birefringence effect. Under such circumstances, the 2f-RT equations 
take a much simplified form which we now state. 

Because f^-', j = 1,2,3,4 are simple (of multiplicity one), expression ([14]) reduces to 

W = 5^iy-D'^(p,p). 

cr 

Consequently (115]) becomes a scalar equation for and the different polarization modes 
decouple: 

(47) iu'W^ - Vx^^" • Vpiy" + Vpfi" ■ VxW^" 

= 27ry"5(^]-(q)-fi-(p))6-(p,q)lr-(q)c^q-2S-(p)ly-(p), Vr 

where 

(48) 6^(p,q) = fi^(p)f]-(q)er(p)<(q)$../,(p-q)e}(p)rf,"(q) 

(49) Wijp) = 7^y"5(^]-(p)-^]-(q))6-(p,q)rfq. 

Note that the Cauchy principal value integral disappears from fl49]) whenever S"^ and W 
commute. 

6.2. Paraxial approximation: spatial anisotropy. Consider now a spatially anisotropic 
spectral density tensor for a medium fluctuating much more slowly in the longitudinal x^ = z 
direction, i.e. replacing $(p — q) in (115]) by 

which, in the limit 6^0, tends to 

(50) S{p-q) J dk^{p^-q^,k). 
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With = ^^(xx, 2, p±,p), the right hand side of eq. (145!) reduces to 

- S-(p^)W-(p^) + W-(p^)S-t(p^^ 
where 6^ = 

'5[aK(P^'^^) = e;'^*(P±,PX'"(q±>p) J '^sifg{p±- (l±,k)dke]:''"{p±,p)dl''^*{q±,p) 



and 
S^(Pj 



5 fi^(p^,p)-fi^(q^,p) 



f]^(p^,p)-n^(q^,p) 



<s^(p±,q±) : Wqj 



Eq. (1451) now takes the paraxial form 

= 2 j 6{n-{p^,p) -n-{q^,p)^e{p^,q^) ■.^W{cl^)dq^ 

(51) - [s-(p^)2n-(p^) + 2n-(px)s-t(p^) 

The longitudinal variable z plays the role of a temporal variable and p is a parameter so 
that (15T|) can be solved as an "initial" value problem given the initial data on z = constant 
and a fixed p if dpQ'^ ^ 0. 

7. Examples 

In this section, we briefly discuss a few media for which the scattering tensor can be 
explicitly computed (see [7] for a more elaborate discussion). 

7.1. Isotropic medium. For the simplest isotropic medium, Kq = diag [eo, Cq, eo, /io; /^O; /^o] • 
There are two nozero eigenvalues: f2'''(p) = co|p|,fi~(p) = — co|p| of multiplicity two. Let 
P = P/ 1 P I s-^d let p\ , p\ be any pair of unit vectors orthogonal to each other and to p so 
that {p\^ pj[, p} form a right-handed coordinate frame. Let {q^, q^, q} be similarly defined. 
The eigenvectors are 



d+'+(p) 



ip± 



ipj 



2P± 
2 P± 



,d-'+(p) 



2 P± 
2 P-L/ 



IP) 



j^P± 
2 P± 



Often, in a scattering atmosphere for instance, /i ~ and consequently 



(52) 



;p, q) 



-<l>,(p-q)pi^q^qi^pl, r, a, z/, C = ±. 



This is the setting for which S. Chandrasekhar originally derived his famous equation of 
transfer In this case, eq. (H5l) is the two-frequency version of Chandrasekhar 's famous 
transfer equation [11[9]. 

As we shall see below, many materials are birefringent and permit two monochromatic 
plane waves with two different linear polarizations and two different velocities to propagate 
in any given direction [3]. This is the birefringence effect. 
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7.2. Chiral media. A chiral medium is a reciprocal, biisotropic medium with the constitu- 
tive matrix 

eol ixi 
-ixl /iol 

where x ^ 1^ is the magneto-electric coefficient. To maintain a positive-definite Kq we 
assume < e/^- We then have 



Kn 



(53) 



PjUjhQ 



Co 


' 


-px 




' zl 


—ml 


1-^2 


px 







inl 





where z = yyUo/eo > is the impedance and k = x^o is the chirality parameter. The four 
non-zero simple eigenvalues and their corresponding eigenvectors are 



-ipi + pI 

—1^1 ■ — 1 " 2 

-z p]_ — iz p]_ 

+ pi 

-z^^p]_ + iz'^p^ 

-tp\ + pi 

z^^p]_ + iz~^p]_ 
ip\_ + pi 

— 1 " 1 • — 1 " 2 

z p]_ — iz p]_ 



Co|p|(l 



K 



,-1. 



co|p|(l - k) 



Co\p\{k - 1) 



Co|p|(-K - 1) 



Note also that = —Q^,Q'^ = —Q^. As \k\ < 1, e^,e2 are the forward propagating modes 
and e^, the backward propagating modes. 

7.3. Birefrigence in anisotropic crystals. The only optically isotropic crystal is the 
cubic crystal. In the system of principal dielectric axes, the permitivity-permeability ten- 
sor of a crystal, which is always a real, symmetric matrix, can be diagonalized as Kq = 
diag[ea;, e^, ez, 1, 1, 1]. One type of anisotropic crystals are the uniaxial crystals for which 
ex = = e± ^ ez = ey (if the distinguished direction, the optic axis, is taken as the z-axis). 
There exist two distinct dispersion relations for the forward modes 



vl + vl 



The backward modes correspond to — fi^, —Vt°. The corresponding wavevector surface con- 
sists of a sphere and an ovaloid, a surface of revolution. The former corresponds to an 
ordinary wave with a velocity independent of the wavevector, the latter an extraordinary 
wave with a velocity depending on the angle between the wavevector and the optic axis [3] . 

Let d°,d^ be the associated eigenvectors. Set Kg = diag[e_L, e_L, ey] and let a'^ solve the 
following symmetric eigenvalue problem: 

(54) - p X (K^) "^p X a" = {n") ^a", a = e, o. 

Then the eigenvectors d*^ can be written as 



(55) 



-p X 
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The same formula applies to the backward modes. Eq. (!54|) has the following solutions 

a' = {-P2,Pi,0y, a° = {pi,p2, 



P3 



from which we deduce that the wave is linearly polarized. 



7.4. Gyrotropic media: magneto-optical effect. For an isotrpic medium [TT] in motion 
or in the presence of a static external magnetic field Hcxt the permittivity tensor Kq is 
no longer symmetrical; it is generally a complex Hermitian matrix. Here we consider the 
simplest such constitutive relation 



(56) 

where g = /Hext, / e 



D = eoE - ig X E, B = H 

is the gyration vector. Equivalently, we can write 



E 



^0 



^ ( eoD + zg X D - -gg^D ) . 

- |g| V eo / 



In this case there are two distinct forward dispersion relations [TT 



n' = c^\p + —g\, 



Co|p-yf 



where Cq = l/^/eQ. Clearly the wave- vector surface consists of two spheres of the same radius 
but different centers. This should be contrasted with the case of chiral media for which the 
wave-vector surface consists of two concentric spheres of different radii. 

The associated eigenvectors d'^, a = 1,2 can be written as in (1551) with a'^ solving (J5^ and 
with 







m 






-m 








m 







Let 
with 



(yfiP^-|-(yf2Pi+5'3P- Wc cau writc the three-dimensional vector a°" as = Pj_+7(tPj_ 

9l -9l- {-lY^/{9i + g^? + Kg^ 



la 



2(9192 - ieo93) 



a 



1,2. 



We see that the wave is in general elliptically polarized or linearly polarized when g is 
orthogonal to the wavevector p and circularly polarized when g is parallel to p. Again, the 
simplicity of the eigenvalues implies that depolarization is absent in the gyrotropic media. 



8. Conclusion 

The main contribution of this work is the derivation of the 2f-RT equations (HSj) . fHTl) for 
the 2f-WD in the weak-disorder regime based on the first Wigner-Moyal equation ([7]). All 
the terms in the equations can be explicitly calculated from the materials properties. 

Let us turn to the second Wigner-Moyal equation (jS]) and briefly discuss its implications. 
By the same multi-scale expansion, the leading order term from eq. is 

2ljW = ipjRjLW + iWLpjRj 
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which, along with (IT^ . then imphes that the wavenumber p should be restricted to the 
surface: 

Denote the area element of the surface by dQ. Hence the mutual coherence in this regime is 
given approximately as 

(u(Xi,ti)u"f(x2,t2)> 



where the coherence matrix W°" = [W^^^] satisfies the 2f-RT equations M5\\ and dQ{p) is the 
area element of the surface uj = Q'^{p). 
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Appendix A. Derivation of Wigner-Moyal equation 
From the Maxwell equations, we have 



2i 



i{27ry 



(27r 



after changing variable and integrating by parts. Using the identity 



Rjd,^ I e-*P'yL(x+^)UiU^(iy 



L(x+^)Ui 



we then obtain 



i—W 



1 



(27r)3 
i 



e-'P'yR,.9.. (L(x+^)Ui lU^rfy 



(57) 

Similarly, 



(5J 



-^-W 



(2vr) 



e-*p'yUiVx ( u5l(x 



2 ■ 



Rjc/y 



£(2vr)3 7 e-^^^^U,UtL(x-|)dyp,R, 
^V.y e-^P^^UiUt2L(x-^)rfyR, 



2(27r) 



By the spectral representation of L we write 

(59) ^ J e-P^^L(x + ^)UiUtdy = J e^^'-L{ci)W{p - ^)rfq 

(60) ^ 



(2vr)= 



e-^p'yUiU^L(x-^)dy 



W(p+ ^)L(q)e*^^^dq. 
2 



Adding or subtracting (1571) and (l58l) with (l59l) - (!60l) we obtain the Wigner-Moyal equations. 
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Appendix B. Derivation of geometrical optics equation 
Consider the following term from V2^'- 

(61) Tr [E"'«"tR.5^^. [L1?„"^D'^'"^]] 

= Tr [W'^^'^YijW^^ [ly^^D"'"^]] + Tr [E'-^'^^Rj-a^Xiy^^D'^'^^] 

= Tr [E"'«^tR,.L5^^^-^D-,"C] + Tr [W^^^^HjlM^^d^p"^''^] + Tr [E^'^^^R^-a^XTy^^D"'"^] 
The first term on the right hand side of (I6T1) can be calculated as 

(62) Tr [E^'^'^tR .L . ly^^^D'^'"^] 

= Tr [e^'-dp^ [e^'«V/R/L] d^W^^^"'"^^] - Tr [e^-^S^^, [e^'^t] p.-RiW^W^^'^'''^] 
= Tr [E^'^'^t^^^^r^^.^^^-^D-.^C] + Tr [e^'^dp^ [e^'^t] [fi- -p.R^L] a,^Ty„^^D'^'<] 

using the eigenvector property and ffTOj) while the last term on the right hand side of fl6T|) is 

(63) Tr [E^'«^tK^.5^. Liy^-^D'^'^] = 5..5,^e^'«tR^.5,Xd^'°iy;,. 

We turn to the middle term on the right hand side of flM]) . We have the following calculation. 

(64) Tr [E^'^'^tRj Liy^^ [9^^ d"'"d"''^"f + d"'°(9^^ d"''^"^] ] 

cr,a,C, 

= e^'«tR^.La,^.d^'"W-;, + 5^ Tr [e^'^^ap, [e^'«V,R,L] lV„^^d'^'"a,,d'^'^t] 

- 5^ Tr [e^'^9,,e^'«tiy^-^p^.R^.Ld-'°9.^d'^'^t] 

cr,a,C, 

= e^^^^n.-Ld^A^^'^Wl, + 5^ Tr [e^-'^Sp, [fi^e^'«t] iy-^d'^'"9,,d'^'^t] 

- ^ Tr [e"''^9p,e"'«tty-^fi-d"'"9^,d"''^t] 

C 

The second part in can be calculated in the same manner. The counterpart of ( !62l) 

yields exactly the same expression as the right hand side of (1^ while the counterparts of 
fl^ and flM|) yield, respectively 

^iTf^d^'^t^^XR^e^'^ 

C 
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and 

a,a,C 

C 

Next let us turn to the first part in PiW: 
(66) Yl [E^'^'^tp^-R, VxL ■ Vp [W^^d'^'^^d'^'^^] ] 

The first term on the right hand side of ( l66i) equals 



while the second term can be calculated as 

a C, a 

(67) + ^ (r^^ - Q'^) W^^d.^e^'^^d'^'^dpA'^'^^e^^'' - J] iy:,a,,e-'«tp .R.Lap^.d^'" 

The last term on the right hand side of (l67j) can be further expressed as 

(68) - Y WLd.,e^'^^d,^ [P.^M^^ + E W^^.^^.e^'^tR^-Ld^'" 

a a 
a a 

+ 5^9,,e^'«tR^.Ld-'°iy:,. 

a 

For the second part in "PiW, the counterpart of (1671) yields exactly the same expression as 
(jnZ]) while the counterpart of (jSZD yields 

5^ e^'^^dpA-'^'^d.^p^W:, + Y W^^d^.n-dpA^^^^^e-'" 

C 

+ ^ (fi^ - ^]'^) ly^^^e^'^t^p^d'^'^d'^'^tS.^e^-'^ + ir^^^d^'^^tLR^.^.^e^''^ - W^^dp^n^d^'^^d^^e^'". 

In the final expression, many of the above terms in PiW — cancel. For instance, all 

the terms involving [Q'^ — Q'^), all the terms involving W^^dp^Q'^, W^j^dp^fl"^ and all the terms 
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involving D.'^W^^, D.'^W^^ cancel. Using the fact e"^'" = Ld'^'" and some algebra we obtain 
(69) Tr [E^'^'^^PiW] - Tr [E'-^'^^PaW] - Vx^^^ • VpW^^ + Vpfi" ■ V^W^^ 

a C C 
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